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Abstract. Weak gravitational lensing has become a powerful probe of large-scale structure 
and cosmological parameters. Precision weak lensing measurements require an understanding 
of the intrinsic alignment of galaxy ellipticities, which can in turn inform models of galaxy 
formation. It is hypothesized that elliptical galaxies align with the background tidal field and 
that this alignment mechanism dominates the correlation between ellipticities on cosmologi- 
cal scales (in the absence of lensing) . We use recent large-scale structure measurements from 
the Sloan Digital Sky Survey to test this picture with several statistics: (1) the correlation 
between ellipticity and galaxy overdensity, Wg^; (2) the intrinsic alignment auto-correlation 
functions; (3) the correlation functions of curl-free, E, and divergence- free, B, modes, the 
latter of which is zero in the linear tidal alignment theory; (4) the alignment correlation func- 
tion, Wg{rp,9), a recently developed statistic that generalizes the galaxy correlation function 
to account for the angle between the galaxy separation vector and the principle axis of el- 
lipticity. We show that recent measurements are largely consistent with the tidal alignment 
model and discuss dependence on galaxy luminosity. In addition, we show that at linear 
order the tidal alignment model predicts that the angular dependence of Wg{rp,6) is simply 
Wg+{rp) cos{29) and that this dependence is consistent with recent measurements. We also 
study how stochastic nonlinear contributions to galaxy ellipticity impact these statistics. We 
find that a significant fraction of the observed LRG ellipticity can be explained by alignment 
with the tidal field on scales > 10 /i~^Mpc. These considerations are relevant to galaxy 
formation and evolution. 
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1 Introduction 

Light travels through the web of cosmic structure on geodesic paths determined by the 
presence of matter. The gravitational deflections due to matter inhomogeneities result in 
distorted observations of distant objects. In special cases, the distortions can be profound, 
producing multiple images and even Einstein rings. More commonly, the distortions are 
subtle and yield only small deviations from the intrinsic shape of the source, a process 
termed "weak lensing." The statistical analysis of these small distortions is becoming an 
important tool in the study of large-scale structure and cosmological parameters [1-4]. Since 
lensing arises from purely gravitational physics, it directly probes the underlying matter 
rather than an observable that only correlates with some fraction of the matter. After 
early successes in detecting the weak lensing signal [5-8], surveys such as COSMOS^ and 
CFHTLS^ have allowed more precise measurements on larger scales. Among other things, 
weak lensing surveys can test General Relativity [10-13], constrain the properties of dark 
energy [14], and measure galaxy bias [15]. Planned surveys such as DES^, LSST^, WFIRST^, 
Euclid^, and Pan-STARRS^ have the potential to further improve the precision of weak 
lensing measurements. 

^ http://cosmos.astro.caltech.edu/; [9] 

^Canada France Hawaii Telescope Legacy Survey, http://www.cfht.hawaii.edu/Science/CFHLS 
^Dark Energy Survey, https://www.darkenergysurvey.org 
*Large Synoptic Survey Telescope, http://www.lsst.org 
^Wide-Field Infrared Survey Telescope, http://wfirst.gsfc.nasa.gov/ 
®http://sci.esa.int/science-e/www/area/index.cfm?fareaid=102 

''Panoramic Survey Telescope & Rapid Response System, http://pan-starrs.ifa.hawaii.edu 
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However, these measurements are far from straightforward, requiring an understand- 
ing of a wide variety of systematic errors (see, e.g., [16] and references therein). The past 
several years have seen a large effort to better understand and correct for instrumental sys- 
tematics. However, even for a perfectly understood instrument, astrophysical uncertainties 
can contribute substantial errors when inferring the matter distribution from weak lensing 
measurements. For instance, the intrinsic orientations and shapes of observed galaxies are a 
source of systematic error. This "intrinsic alignment" (lA) of galaxies is particularly prob- 
lematic since it can be significant and would bias even an ideal measurement. Weak lensing 
statistics involve averaging across pairs of observed galaxy ellipticities 7°'^'* , which consist of 
both the intrinsic ellipticity (I) of the galaxy and the gravitational lensing shear distortion 
(G): 7°*^*^ = 7^ -|- 7*^. For the i-component of ellipticity, the observed ellipticity correlation 
function is then 

^^obs^obs^ = (7f 7f ) + {ih?) + {ifiD + {ihD- (1-1) 

The desired signal in weak lensing studies is the first term (GG). If intrinsic alignments 
are random, the GI and H terms average to zero. However, the coherent influence of large- 
scale structure on galaxy ellipticity contaminates the lensing signal because intrinsic shapes 
acquire a non-zero average correlation. Since the weak lensing signal is small (7*^ is roughly 
1% of 7^ for a typical galaxy [17]), even small correlations can lead to appreciable intrinsic 
alignment contributions. It has long been known that the II correlations from galaxies in 
close proximity (where intrinsic ellipticity correlations are strongest) could be a significant 
contaminant [18-22]. Fortunately, the II term can be easily reduced by either down- weighting 
or excluding nearby pairs [23-26]. It was later realized by [17] that the GI term can also 
be a significant contaminant, introducing a correlation in the ellipticities of objects that are 
along the same line-of-sight but separated by a large spatial distance. A foreground lensing 
potential affects the intrinsic ellipticities of nearby objects as well as the observed ellipticites 
of background objects via lensing. Observations have confirmed the presence of both of these 
intrinsic effects [16, 27-30]. 

It is critical to understand lA for high-precision weak lensing experiments. In [27], it was 
shown that GI contamination at the ~ 10% level is possible in a typical weak lensing survey. 
The potential degradation of cosmological parameter measurements by lA contamination is 
significant. For instance, lA can bias cosmic shear measurements of a% at the current level of 
uncertainty (ag ~ 0.8 it 0.07) for a CFHLTS-like survey, which is larger by ~ 3 than the best 
current constraints [31]. Similarly, uncertainty in the amplitude of intrinsic alignments can 
impart a significant bias in cosmological parameter measurements, even when a particular 
model is assumed in order to subtract the alignment signal [32]. However, these alignment 
effects are not just a contaminant - they also provide a probe of large-scale structure and 
galaxy formation. 

Several models of lA with varying levels of complexity have been proposed (e.g. [20, 
31, 33, 34]). These models belong to two general classes: alignment of the galaxies with 
background tidal field or torquing of the galactic angular momentum vector by the tidal field. 
Any analytic predictions of the orientation and ellipticity of a galaxy residing in a background 
tidal field relies on assumptions relating the orientation of dark matter halos with that of 
the resident galaxies, and nonlinear scales are particularly difficult to model analytically. 
Intrinsic alignment of dark matter halos has also been studied using A^-body simulations 
[28, 35-38], which have shown that ellipticities can differ significantly in both shape and 
orientation when measured in the inner and outer regions of halos. A better understanding 
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of the relationship between halo and galaxy ellipticities may require simulations with baryon 
physics and an exceptional dynamic range (see, e.g., [39]). 

This paper tests analytic models and different statistical measures. In particular, we 
focus on the linear tidal alignment (LA) model [17, 20], which posits that the intrinsic 
ellipticity of a galaxy is a linear function of the tidal field. This model should dominate 
on large scales for elliptical galaxies. Recent work by [32] has shown that the model is 
consistent with measurements of GI correlation. We expand this comparison and consider 
possible extensions to the model. 

This paper is organized as follows. Section 2 discusses intrinsic alignment models, 
especially the linear tidal alignment model. In Section 3, we summarize recent measurements 
of intrinsic alignment. We then review several statistics relevant for GI and II correlations 
and calculate these statistics in the LA model. One of these statistics is the alignment 
correlation function, Wg{rp,9), recently proposed in [28]. By adding an angular dimension 
to the correlation function, this statistic can in principle contain additional information on 
the relationship between galaxy clustering and alignment, a prospect which we examine. 
We compare the model predictions to recent measurements and determine the consistency 
and strength of linear tidal alignment. We also propose a potential signature of nonlinear 
alignment contributions and nonlinearities in the density field. In Section 4, we consider 
the effects of a stochastic contribution to galaxy ellipticity that does not correlate between 
galaxies. We study two models for this stochastic component and discuss the impact on the 
measured alignment statistics. Section 5 summarizes our conclusions and provides discussion 
in the context of galaxy formation. We also include an appendix with the details of some 
calculations referenced in the text. Throughout this work, we assume a flat ACDM cosmology 
with Qrn = 0.25, Qb = 0.04, as = 0.8, = 1, and Hq = 70 km s"^ Mpc'^. 

2 The linear tidal alignment model 

The collapse of overdense regions into dark matter halos and galaxies occurs preferentially 
along the stretching axis of a background tidal field, and galaxy intrinsic ellipticity should 
maintain some memory of this asymmetry at the time of formation [20]. In particular, 
elliptical galaxies are supported primarily through velocity dispersion rather than rotation 
and are thus more likely to align with the surrounding halo and background tidal field. It 
is physically reasonable for galaxy orientation to correlate with the principal axis of the 
gravitational tidal field [17]. The LA model of [20] relates the intrinsic ellipticity* of an 
elliptical galaxy to a linear function of the tidal field: 

7(+,x) = -^(^' - V2,2V,V,)5[^p], (2.1) 

where Ci parameterizes the strength of the alignment, with sign convention such that positive 
Ci corresponds to preferential galaxy alignment along the stretching axis of the tidal field. ^ 

^The quantity used here is actually the intrinsic shear, which differs from ellipticity by a factor 1/2R, 
described after eq. 2.2. To avoid confusion with lensing shear, we refer only to intrinsic ellipticity. 

^We define Ci to capture the full magnitude of the LA effect. In several previous studies (e.g. [32]), Ci 
was specified using a standard but somewhat arbitrary normalization calculated from ellipticity variance [17], 
and an additional dimensionless constant parameterized the strength of LA with respect to this reference 
value. As we discuss in section 4, stochastic contributions can affect large-scale correlations differently than 
ellipticity variance. Since the LA model is most applicable on large scales, we choose a convention for Ci that 
relates it directly to the magnitude of large-scale correlations. 
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Positive Ci thus yields an anti-correlation between the intrinsic alignment of a foreground 
object and the gravitational shear of a background object, ^'p is the gravitational potential, 
and S is a filter that smooths fluctuations on halo scales. In [17], S is chosen to be a top- 
hat in Fourier space with a maximum wavevector of 1 /iMpc~^. When compared with no 
smoothing, we find that this choice has a negligible effect on scales of interest (> 10 /i~^Mpc), 
and we thus effectively ignore S. Up to derivatives, eq. (2.1) is the unique function of 
that is local, linear, and quadrupole symmetric. Since higher-derivative terms should be 
negligible on large scales, the LA model is unique up to the normalization Ci. The x- and 
y-axes in eq. (2.1) are on the plane of the sky, and ellipticity is decomposed with respect to 
this coordinate system: 
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where (j) is the position angle measured from the x-axis, and h/a is the axis ratio. R is the 
shear responsivity, which captures the average response of measured ellipticity to a small 
shear [40]. 

The density- weighted intrinsic ellipticity is defined as 7^ = (1 -|- Sg)^^ , for galaxy over- 
density 5g = bg5, where 5 is the matter overdensity, and bg is the galaxy linear bias factor. 
In the linear regime, the gravitational potential at redshift Zp is related to the density field 
via the Poisson equation (valid on sub-horizon scales, k » cHq): 



^p(fc) = -AnGpmfiil + zp)k 2(5(fe, zp), 



(2.3) 



where pm.o is the present mean matter density. There is some ambiguity to the appropriate 
redshift, zp, to evaluate the potential. In previous work [17, 32], zp was chosen to be 
during matter domination, when most of the stars in elliptical galaxies formed. However, 
it is also plausible that recent accretion significantly impacts the alignment, in which case 
the gravitational potential should be evaluated at roughly the observed redshift. In linear 
theory, the only difference between the potential at different times is the overall amplitude, 
which can be absorbed into Ci. Thus, the inferred strength of lA will depend on zp. For 
consistency with previous work, we evaluate ^p during matter domination. 

A product of density fields in configuration space becomes a convolution in Fourier 
space, and thus the linear model predicts: 
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where k2 = k — ki, 6^^\k) denotes the 3-dimensional Dirac delta function, and D{z) is the 
growth factor, normalized so that (1 + z)D(z) = 1 during matter domination. 

The LA model breaks down on nonlinear scales. As in recent work [31, 32], we estimate 
these effects by using the Halofit nonlinear density power spectrum [41]. We note, however, 
that simply applying a nonlinear power spectrum to the LA model is not a full nonlinear 
theory. Moreover, the resulting nonlinear corrections depend on the choice of zp, since the 
growth function does not fully capture the evolution of the nonlinear density field. For 
instance, choosing zp at the observed redshift (which provides the maximum nonlinear cor- 
rection) rather than during matter domination can affect the predicted alignment amplitude 
by ~ 20% for GI correlations at scales of 5/i~^Mpc. 
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NL corrections improve fit below 10 /i^^Mpc 




0.123 ±0.014 
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NL corrections improve fit below 10 /i~^Mpc 
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WE 


Use 10++ fit 
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0.02 
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0.68 


0.64 


LA prediction is wb = 




Cl Pcrit 








Wg+ 


0.71 ±0.02 


1.8 


0.12 


Calculated without weighting by 70 


W+ + 


0.74 ±0.07 


0.24 


0.91 


Calculated without weighting by 70 


Wg{rp,9) 


0.16 - 1.55 






Luminosity dependent - see section 3.4 



Table 1. Summary of tests of LA model, including measured model parameter and reduced ^oi 
the fit. All measurements use the SDSS LRG catalog except for Wg{rp,d), which uses the SDSS main 
sample. The upper section contains statistics that are weighted by the ellipticity magnitude while the 
lower section contains unweighted statistics. 

On large scales, the LA mechanism should dominate ellipticity correlations, which will 
scale linearly with the matter power spectrum, Ps{k). Spiral galaxies are supported by an- 
gular momentum, and thus a distinct alignment mechanism, based on the tidal torquing 
theory of protogalaxies, may be relevant. Models based on tidal torquing can be categorized 
as "quadratic alignment models," since the tidal field enters quadratically at lowest order 
rather than linearly [20, 34, 42], suppressing large-scale correlations because 6 <^ 1. Predic- 
tions of intrinsic alignment effects from quadratic models are qualitatively different from the 
linear model. For example, quadratic models predict a divergence-free (B-mode) component 
to the ellipticity at leading order [43] but a vanishing lowest-order correlation between matter 
density and ellipticity. Nonlinearities in the density field could potentially allow quadratic 
alignment effects to contribute at linear order in Ps{k) [34]. Recent observations [27, 28, 44] 
have split galaxies by color into "red" and "blue" sub-samples, finding qualitative differences 
in intrinsic alignment, suggesting the possibility of different alignment mechanisms. Blue 
samples exhibit weaker intrinsic alignment on large scales, supporting the theory that LA 
effects are less prominent in spirals. 

3 Measuring intrinsic alignment 

There are numerous probes of galaxy intrinsic alignment. We consider several alignment 
statistics in real space (for both GI and II correlations) and compare measurements with LA 
model predictions. Table 1 provides a summary of these statistics, which are described in 
the following subsections. 

3.1 Galaxy samples 

To test the predictions of the linear alignment model, we compare with existing measure- 
ments of intrinsic alignment statistics. Catalogs from large, deep surveys have allowed recent 
measurements of these correlations with better precision and at larger separations than was 
previously possible. In [29, 30], the authors used the catalog of Luminous Red Galaxies 
(LRGs) from DR6 of the Sloan Digital Sky Survey (SDSS; [45]) to measure both II and GI 
correlations. LRGs, which are among the most luminous elliptical galaxies, are expected to be 
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particularly well described by the LA model. The LRG sample used in their analysis contains 
83,773 objects with spectroscopic redshifts in the range 0.16 < z < 0.47 and mean redshift 
of z = 0.32. Note that the LRG shapes used in [29, 30] are measured without correcting 
for the point-spread function (PSF). We find that 70 is larger by 10-30% when measured 
without PSF-correction (R. Mandelbaum, private communication). This offset would affect 
the amplitude of any lA statistic that weights by ellipticity magnitude and thus the inferred 
value of Ci. 

We also consider the measurements made in [28] using SDSS DR6 main sample galaxies 
from the New York University Value Added Galaxy Catalog [48]. This sample includes 
430,164 galaxies with spectroscopic redshifts 0.01 < z < 0.40. For more information on the 
redshift distribution of this sample, divided by luminosity, see [47]. 

Since we wish to test the LA model through comparison with observations of the el- 
lipticity statistics, we face the opposite problem as cosmic shear experiments: weak lensing 
introduces a contaminating signal into the observations. However, the lensing signal is neg- 
ligible for the measurements employed here, which have low mean redshifts and where only 
correlations between galaxies within a projection volume of ~ 100 Mpc are included. 

3.2 Intrinsic alignment statistics 

Intrinsic alignment contributes both a cross-correlation (GI) term and an auto-correlation 
(II) term (see [16, 17, 27]). We write the intrinsic ellipticity cross- and auto-correlation 
functions as: 



where i = {+, x}, and these two components are measured with respect to the separation 
vector on the sky. Note that by symmetry, (f) = 0. Since it is the background matter field 
that is responsible for gravitational lensing, S^g^ provides a measure of the GI correlation, 
while probes the II correlations. 

The separation vector r can be separated into a component on the sky (rp) and along 
the line-of-sight (11). The projected correlation function of quantity X can be written in 
terms of the 3D correlation function as 



where 2nniax is the depth of the projected volume. 

It is straightforward to calculate these correlation functions in the LA model. Unless 
otherwise specified, we include terms to quadratic order in 6 in the following calculations. In 
the LA model. 



^g,ir) = {5g{x)^Ux + r)) 
Cn{r) = {li{x)^l{x + r)) 



(3.1) 




(3.2) 



(5,(0)7i(rp,n)) 




(3.3) 
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We have defined the x-axis to be the separation axis on the sky and the z-axis to be along 
the hne-of-sight. n is the magnitude of the wavevector on the sky, while is the magnitude 
along the line-of-sight (fc^ = + k1). The projection of ^^_J^ along the line-of-sight is 

Wg^{rp) = ~| ^^^^'° j dKdkz^^Ps{k,z)sin{kzU^s^)J2{K,rp). (3.4) 
Similarly, the projected auto-correlation statistics are given by 

w^f_fl,_,xx)(^p) = ^ ( D{z) ) J '^'^'^^^■p^-f'5(^'^)sin(A;2nmax) [Mt^rp] ± J4(Krp)] . 

(3.5) 



In figures 1-2, we compare the linear tidal alignment model predictions with measure- 
ments from SDSS LRGs [29, 30], which project these statistics over 

Hmax = 80 /j-^Mpc. Er- 
rors were calculated through jackknife resampling. Redshift-space distortions are included in 
the model predictions, as described in Appendix C. We measure the galaxy bias bg by fitting 
to the measured projected galaxy correlation function, Wg{rp), and find 6^^*^ = 2.12 =b 0.04. 
The II correlations presented in [29] are functions of 3-dimensional separation rather than 
projected separation. To compare with the LA prediction, we have assumed that the statis- 
tics are isotropic and performed a projection along the line-of-sight following eq. (3.2). To 
calculate the errors for these projections, we use the mean and variance of the 3-dimensional 
measurements and perform the projection on 1000 random realizations. In addition, the esti- 
mator used to calculate the II statistics in [29] results in an additional factor of (1 -|- Cg('"))~^) 
where £,g{r) is the 3-dimensional galaxy correlation function. We have removed this factor 
before performing the projection. For reference, the open circles in the left panel of figure 2 
show the measurements without this correction. 

To calculate the magnitude of Ci, we fit to the data at large separations (rp ~ 10- 
140/i~^Mpc), where we expect the linear model to hold with minimal contamination from 
other alignment effects. Below this scale, nonlinear effects become appreciable, and the LA 
model is no longer valid. The LA theory agrees quite well with measurements at separations 
above ^ 20 h~^Mpc for Wg+ (figure 1) and above 10 /i~^Mpc for w++ (figure 2). Including 
nonlinear corrections with Halofit improves agreement on smaller scales. Agreement is weaker 
but still reasonable for the Wxx statistic. 

Table 1 shows the fit results as well as reduced to indicate the goodness of fit. We 
quote values for the dimensionless quantity Cipcrit, where pm,o = ^mPcrit- We fit separately 
to the II and GI statistics, both with and without ellipticity weighting. Calculating these 
statistics without weighting by ellipticity is equivalent to setting the ellipticity magnitude, 
7o, of each galaxy equal to 1. We define Ci, the non- weighted analog of Ci. The ratio of 
weighted to non- weighted statistics is ~ 0.17, although as discussed below, some care must 
be taken to distinguish between this value and the mean ellipticity of the LRG sample. 

We find that the w++ and Wg+ statistics are well-described on large scales by the LA 
model with consistent amplitudes. The value of Ci measured here is consistent with that 
found in [32], who examined 'Wg-^- for a variety of data sets and found Cipcrit ~ 0.13 it 0.02 
for a similar SDSS LRG sample (medium luminosity bin).^*^ 

^"Note that [32] use PSF-correctcd galaxy shapes. 
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Figure 1. Measurements of [30] and LA model prediction for The black dashed line is calculated 
using the linear theory Ps{k), and the red solid line uses the Halofit model. 




Figure 2. Measurements of [29] and model predictions for w++ {left panel) and Wxx {right panel). 
The measurements have been projected along the line-of-sight. Open circles, indicating the original 
measurements without the (1 + ^g{r)) correction, are only shown for w++ and on small scales where 
there is an appreciable difference. For clarity, these points have a small horizontal offset. Line 
convention is the same as in figure 1. A linear y-axis is used for Wxx- The normalization of the LA 
prediction for both statistics is set from the fit to 



3.3 Autocorrelation E- and Z?-modes 

The Wxx and w++ statistics can be written in terms of curl-free {E) and divergence-free {B) 
modes. Lensing by matter produces only E'-modes, making such a decomposition a useful 
diagnostic in studying the effects of intrinsic alignment and other systematics [49] . As shown 
below, only i?-modes are produced in the LA model, and thus i?-modes could indicate the 
presence of separate alignment mechanisms [43]. 

Following [50] , we can express the E- and S-components of the auto-correlation func- 
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Figures. Left panel: i?-mode auto-correlation statistic. Line convention and LA model normalization 
are the same as in figure 2. Inset shows more detail above = 10 ft,~"'^Mpc. Right panel: _B-modc 
auto-correlation statistic is compared with the LA prediction of zero. The observations are consistent 
with the prediction above 10 /i~^Mpc. 

tions in real space in terms of the linear combinations w±{rp) = w^+{rp) ± Wxxifp): 

w+{rp) ±w'{rp) 
W{E,B){rp) = ^ —, (3.6) 

where w'{rp) is a non-local function of W-{rp): 

w [rp) = w-(rp) +4 dr 12rp / dr —. (3.7) 

J Tp ^ J Vp ^ 

Prom the integral properties and recursion relations of Bessel functions (see Appendix B), 
w'{rp) = w+{rp) in the LA model, which allows us to write the E/B decomposition: 

'i"E^{rp) = J dKdkz-^^Ps{k)sin{kJlmax)Jo{K-rp), 

w'^B^rp) = 0. (3.8) 

Projecting the auto-correlation measurements of SDSS LRGs and using eqs. (3.6) and 
(3.7), we have calculated the observed E- and S-mode signals. Errors are calculated as with 
wj^j^ and Wxx- Figure 3 shows the LA predictions and measurements. The observed wb{tp) 
is consistent with the LA prediction of zero on scales above 10 /i~^Mpc. Below this scale, 
nonlinear terms become important, and it is not expected that wb would remain negligible. 

3.4 The alignment correlation function 

Faltenbacher et al. [28] introduced the alignment correlation function, ■Wg{rp,9), which de- 
scribes the dependence of clustering on both projected separation, Vp, and the galaxy orien- 
tation angle, 9, measured from the axis of separation. It is related to the previously defined 
correlation statistics via the relations 

Wg{rp) = — / d9wg{rp,0); Wgj^ = — I d9 cos{29)wg{rp,9), (3.9) 
71" Jo Jo 



-9- 



where Wg+ = {cos{2(j))5g) is the Wg+ statistic when not weighted by elhpticity. Because of its 
angular dependence, w^{rp,9) can in principle provide more information than Wg+. Indeed, 
eq. (3.9) shows that Wgj^ is the dipole moment of Wg{rp,9) and thus lacks any information 
from higher angular moments. 

We now derive Wg{rp,9) in the LA model. An arbitrary periodic function of 9 can be 
written as a sum of cos{n9) and sm{n9) terms. However, two symmetries exist for Wg{rp,9) 
which greatly restrict possible terms. First, by parity, Wg{rp,9) = Wg{rp, —9), ruling out all 
sin(n^) terms. Second, since the ellipticities are invariant under rotations by vr, Wg{rp,9) = 
Wg{rp, 9-\-tt), which allows only cos{n9) terms with n even. In its most general form, therefore, 

Wgirp, 9) = Wg{rp) + ^ a„(rp) cos(n6'). (3.10) 



In Appendix A, we demonstrate that the only possible angular dependence when density 
and elhpticity fields are Gaussian is the cos(2^)-term. Under these conditions, which are met 
in the LA model, and applying eq. (3.9), we find: 

Wg{rp, 9) = Wg{rp) + 2wg+{rp) cos(20). (3.11) 

Thus, Wg{rp,9) contains the same information as Wg{rp) and Wg-^{rp) in the LA model. 

Figure 4 compares the LA model prediction of the alignment correlation function with 
the SDSS main sample measurements of [28]. These measurements divide galaxies by lumi- 
nosity and type. Errors represent the variance between 10 realizations in which galaxy orien- 
tations have been randomly shuffled, as described in [28]. We consider red galaxies in the fol- 
lowing four magnitude bins: {-20 < Mr < -19; -21 < Mr < -20; -22 < M,. < -21; -23 < 
Mj. < —22}. Consistent with the measurement, a projection length of Ilmax = 

40/i~^Mpc 

is used to calculate the LA model prediction. Galaxy bias for each luminosity bin is calcu- 
lated by fitting to the observed correlation function (above 10 /i~^Mpc) and assuming that 
the total sample has bg = 1.07 it 0.01 [51]. The best-fit values of LA strength, in order of 
increasing luminosity, are CiPcrit = {0.16 ± 0.01; 0.17 ± 0.02; 0.36 ± 0.03; 1.55 ± 0.16}. We 
speculate on what the trend of stronger alignment for more luminous objects may imply 
about the dynamics of galaxy alignment in section 5. The solid curves in figure 4 include 
redshift-space distortions, while the dashed curves do not. Redshift-space distortions alter 
the ratio Wg{rp,9) /wg{rp) by over 50% on large scales. The galaxy correlation function is 
more enhanced than the angular term because it involves an integral over Jo^KVp) rather than 
J2{Krp) and thus receives greater contributions where enhancement is large (see appendix C). 

3.4.1 Weighting by galaxy ellipticity 

Recent work [30] has detected a small correlation between 9 and 70 that affects the measured 
ratio between ellipticity-weighted and unweighted GI statistics on all scales. In the absence of 
such a correlation, the ratio between ellipticity-weighted and unweighted statistics is simply 
the sample mean ellipticity (70)- This correlation can arise, for instance, if galaxies pointing 
towards over-dense regions have larger average ellipticities. This ratio depends on correlations 
of both 70 and 6g with 9. For Gaussian density and ellipticity fields, we find that the effect of 
this correlation is to modulate the ratio at all separations by a factor A/tt (see Appendix A): 

= — = z(7o) = (3.12) 



%+ Ci vr 
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Figure 4. Top panel: The best-fit LA model prediction for Wg{rp,d)/wg{rp) compared with the 
measurements of [28] for 4 luminosity bins. The Halofit Ps{k) is used. Solid curves include redshift- 
space distortions, while dashed curves do not. The measurements are divided into 3 angular bins. 
Bottom panel: Deviations from the prediction of cos(20) angular dependence. See section 4.2 for more 
information. For clarity, small horizontal offsets have been introduced to circle and diamond points. 



where 7rms = y (To)- "^^^ ~ 15 it 10% excess of this ratio found at all separations in [30] 
is roughly consistent with the A/tt — 1 = 27% excess predicted by the LA model. However, 
alignment stochasticity may be able to suppress this ratio (section 4). 

3.4.2 Higher-order angular dependence 

Detection of cos{n9) angular dependence in Wg{rp, 9) for n > 4 at a given scale would indicate 
contributions from nonlinearities (which introduce non-Gaussianity in the density field) or 
other alignment mechanisms. 

Detecting such dependence is challenging, further complicated by the fact that small- 
scale alignment processes suppress these terms (see section 4). We can fit observed or simu- 
lated data to a model which includes a higher-n angular term: 

w{rp, 9) = w{rp) + 2wgj^{rp) cos(2^) + W4^{rp) cos(46'). (3.13) 

In the case of the binned alignment correlation function presented in [28], the low angular 
resolution makes a reliable fit challenging. Nevertheless, one can look for a potential cos(40) 
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term by comparing observations with the LA prediction. In the absence of higher-n terms, 



2u;g(rp,7r/6 < 6* < 7r/3) = Wg{rp,0 <e < tt/Q) + Wg{rp,7r/3 <e < 7r/2) = 2wg{rp). (3.14) 

In the bottom panel of figure 4, we show that measured deviations from this prediction are 
largely consistent with zero, although the errors are large. Future observations with improved 
precision and better resolution in alignment angle should allow for detection of (or constraints 
on) higher-n terms arising from processes outside the LA model. This analysis can also be 
applied to halo orientations from simulations. 



4 Additional contributions to intrinsic ellipticity 

The LA model as described in section 2 assumes that tidal alignment occurs down to some 
minimum scale, determined by the smoothing filter S, and processes occurring on smaller 
scales are ignored. We now consider extending the model by including alignment stochas- 
ticity due to the small-scale processes. For example, this stochasticity could arise from the 
small-scale tidal field or internal processes within the galaxy that generate ellipticity. Such 
stochasticity need not be small and will result in galaxy ellipticities that deviate from the 
predictions of any large-scale model in an uncorrelated fashion. The presence of such stochas- 
ticity is physically motivated, and [28, 29] have shown that misalignment between halo and 
galaxy ellipticities is necessary to match the observed alignment amplitude. 

We first consider a model in which stochasticity results in a Gaussian-distributed mis- 
alignment angle 6m with width a. This model has been previously used [29, 30] to introduce 
misalignment between the orientations of simulated halos and the galaxies that populate 
these halos. The full misalignment between galaxies and the background tidal field is the 
combination of misalignments between (1) the tidal field and halo and (2) halo and galaxy. 
Different physical processes may contribute to each. We define the probability distribution 
of the misalignment angle, 6m, with width a: 



f{6m)d6r, 



1 



27ra 



exp 



d6r^ 



(4.1) 



We now calculate the effect of this misalignment on the alignment statistics. The true 
alignment correlation function is the convolution of the original Wg{rp,6) and the misalign- 



w 



true 



{rp,6) = Wg{n 



f{6). This convolution provides the effect of 



ment distribution: 

stochasticity on which is simply the dipole moment of Wg(rp,6). Angular dependence 
in Wg{rp,6) appears as cos{n6) terms, so the relevant convolution is 



/ 



d6m cos{n{6 - 6m)) 



2-Ka 



exp 



cos(n^) exp 



1 2 2 
--n a 

2 



(4.2) 



For each cos{n6) term, this model leads to a stochastic suppression factor of S'n, stoc = 
exp [— in^cr^] . Stochastic angular misalignment suppresses the lA signal and does so particu- 
larly strongly for higher-order effects (i.e. for larger values of n). Physically, this exponential 
suppression comes from periodicity under rotation: a stochastic rotation of ~ vr/n erases the 
signal. Since Wgj^ comes from the n = 2 term of the alignment correlation function, stochas- 
ticity provides a suppression of 52,3100 = exp [— 2cr^] . Similarly, and Wxx come from the 
n=2 term of an analogous alignment correlation function with two independent angles which 
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both contribute in the convolution. Xhus, they will be suppressed by 5*2 siq,^ — exp [-4^2]. In 
the previous sections, the effects of stochasticity are included in the values of Ci determined 
by fitting to measurement. The magnitude of Ci thus reflects both the degree to which 
galaxies align with the tidal field and the amount which they are stretched along it. 

The analysis of [29] measured the misalignment between halos in an A^-body simulation 
and the resident galaxies required to match the observed alignment amplitude, finding a ~ 
35° = 0.61 rad. This value ignores misalignment between the background tidal field and the 
halos and thus provides a lower limit to the total misalignment of galaxies. Using this value, 
'S'2,stoc ~ 0.5, and the suppression of an n = 4 term is sslO times greater. 

The angular-misalignment model discussed above does not include any stochastic con- 
tribution to the magnitude of ellipticity and is thus more applicable to the unweighted statis- 
tics, which do not depend on this magnitude. To include more general stochastic effects, we 
consider a second model that assumes a Gaussian distributed scatter in 7+ and 7x. For a 
distribution width of A^, this yields stochastic contributions of 7+^ and jxs with 



/(7+s,7xs)c?7+sd7> 



1 



7rA2 



exp 



(4.3) 



Note that this model is equivalent to adding Gaussian scatter in the perpendicular com- 
ponents of observed galaxy shape. The effect of stochasticity in this model is to increase the 
total value of galaxy ellipticity: 7^^^^, ^^^^^ = 7^j^g -|- A^. Since it adds only a non-correlated 
component to ellipticity, this model does not affect the magnitude of Wg^ relative to the LA 
model with no stochasticity: {6 {'j^^ + 7+°^)) = (57+^)- However, by increasing 7rms,totaii 
it suppresses Wg^ and thus the cos(20)-term of Wg{rp,9). The stochastic contributions to 
7+ and 7x result in both a change in the observed ellipticity magnitude and alignment an- 
gle. The suppression of Wg^ can be seen as coming from the angular misalignment, which is 
qualitatively similar to the angular misalignment model discussed above. Both models thus 
result in suppressed angular information, particularly for higher-n terms. 

In addition to affecting the alignment correlation statistics, astrophysical stochasticity 
alters the observed value of 7rms- In principle, one can measure the magnitude of stochas- 
ticity by comparing measurements of 7rms with the LA theory prediction in the absence of 
stochasticity: 7rms la ~ — 0)- The true 7rms will exceed the stochasticity-free pre- 

diction. Although stochasticity in the angular-misalignment model makes no contribution to 
7rms; it suppresses the measured value of Ci by S'2,stoc, leading to an inferred 7rms less than 
the true value. However, it is challenging to draw conclusions from this comparison because 
the integral for .^Jjl^(O) is highly dependent on the smallest scales for which the LA model is 
assumed to be valid (i.e. the choice of smoothing filter). This dependence is illustrated in 
figure 5. For instance, smoothing on a typical halo scale corresponds to a maximum wavevec- 
tor fcmax ~ 27r/l /iMpc~^, while smoothing on nonlinear scales corresponds to /Cmax ~ 27r/10 
h Mpc~^. Due to these complications, we have chosen to leave 7rms as a parameter to be mea- 
sured independently. The value of Ci can then be determined from the alignment statistics 
rather than from 7rms; as has been sometimes done in the absence of alignment measurements 
(e.g. [17, 20]). Nevertheless, as seen in figure 5, alignment with the large-scale tidal field can 
account for a large fraction of the total observed LRG ellipticity. Note that for k^ax ^ 1 ~ 10 
/iMpc~^, the LA model will start to break down and the predictions shown here will be of 
limited validity. 

Finally, we consider whether stochasticity suppresses the ratio Wg+/wg+, discussed in 
§3.4. In general, stochasticity can suppress this ratio, which depends on correlations of both 
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Figure 5. The LA model prediction for 7i.nis plotted as a function of the cut-off value fcmax- Solid 
lines are calculated using the Halofit Ps{k), while the dashed lines are calculated with the linear 
Ps{k). The prediction is normalized assuming the best-fit LRG value for Ci. Black lines assume a 
level of stochasticity consistent with [29], while red lines assume that no stochasticity is present. For 
reference, the observed LRG value of 7rms ~ 0.17 is plotted (horizontal blue dotted line). 

7o and 5g with 9. However, in the two models considered in this section, the ratio is unaltered 
because any decrease in the correlation between 70 and 6 is exactly offset by a decrease in 
the correlation between 5g and 6. 

5 Summary and conclusions 

This work examined intrinsic galaxy alignment in the context of the tidal alignment model. 
Using both II and GI correlation statistics in configuration space, we showed that the model is 
consistent with measurements on scales above > 10 /i~^Mpc. Including nonlinear corrections 
in Ps{k) improves the fit at smaller separations, although we note that this is not a rigorous 
nonlinear model. We measured best-fit values of the model strength parameter Ci using 
several statistics and two data sets, finding results that are consistent both internally and 
with previous measurements. We examined lA statistics both unweighted and weighted by 
ellipticity magnitude, 70. The ratio of these statistics (e.g. Wg+/wg+) depends not only on 
the mean value of ellipticity for the observed galaxy sample, (70), but also on correlations 
between both 70 and 6g and the position angle, 6. The LA model makes predictions for these 
correlations which are qualitatively similar to what was observed in [30]. 

The LA model predicts zero i?-mode auto-correlation signal at lowest order. We found 
that LRG measurements are consistent with this prediction. More precise measurements of 
the i?-mode signal should provide a valuable way to determine the importance of higher-order 
terms in 5 in the LA model, other alignment mechanisms, and systematic errors. 

The recently proposed alignment correlation function Wg{rp,9) is of particular inter- 
est as its angular dependence can contain additional information on galaxy alignment and 
large-scale structure. However, we showed that for Gaussian density and ellipticity fields, a 
condition satisfied in the LA model, Wg{rp,0) has angular dependence cos(20) and is com- 
pletely described by Wg{rp) and Wg-^-{rp). This angular dependence is consistent with the 
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recent measurement of [28] . We found that departures from these conditions allow additional 
angular terms to contribute, although they are suppressed by alignment stochasticity. The 
next order angular term is suppressed by an order of magnitude for the level of misalignment 
that has been measured [29, 30]. 

Finally, we considered how to include stochastic contributions to a model of intrinsic 
alignment. Astrophysical sources of stochasticity can affect both the magnitude and direction 
of galaxy ellipticity, and previous measurements [28-30] indicate that such stochasticity is 
needed to account for differences between the strength of alignment measured in simulated 
halos as opposed to observed galaxies. This stochasticity is compounded by uncertainties in 
measuring (70): correcting for the PSF can alter the observed value by ^ 10-30%. Due to 
these complications, it is not recommended to normalize any lA model using local statistics 
such as ellipticity mean or variance as some studies have done. Nevertheless, the LA model 
is able to account for a significant fraction of observed LRG ellipticity. We considered two 
models for stochasticity, angular misalignment and ellipticity scatter, and found that these 
models have similar observational signatures and are consistent with current measurements 
of intrinsic alignment amplitude. Since much of this stochasticity should come from physical 
processes such as galaxy mergers and gas dynamics on nonlinear scales, it is likely that 
the stochastic suppression will be a function of time, and thus care must be taken when 
attempting to separate the intrinsic alignment and weak lensing signals based on their redshift 
dependence (e.g. [20]). 

Previous work (e.g. [28, 32]) has found a significant luminosity dependence in the 
strength of galaxy alignment. The tendency of stochasticity to vary with time will affect the 
relative amplitudes of alignment correlation statistics between galaxy samples that formed 
at different redshifts. One possible explanation for the trend of alignment increasing with 
galaxy luminosity is that more luminous objects have formed more recently and have had less 
time to misalign from the tidal axis along which they formed. This picture is consistent with 
a scenario in which halo ellipticities are determined primarily during times of collapse and 
accretion. For instance, one plausible mechanism for forming triaxial halos aligned with the 
tidal axis is the radial orbit instability (ROI; see, e.g., [52, 53]). The ROI occurs when par- 
ticles with radial orbits are torqued by a small asymmetry in an otherwise spherical system, 
causing them to align with and enhance the asymmetry. The ROI is potentially responsible 
for the overall triaxial shape as well as the nearly universal radial density profile of halos 
[52-54]. We speculate that the ROI (or some other mechanism) acts during the accretion 
of matter to align the halo and the galaxy it contains with the tidal field. Such a mecha- 
nism may allow the accretion of a small fraction of total halo mass to affect the shape and 
orientation of the entire halo. 

However, this alignment mechanism is likely a transient phenomenon. After accretion 
stops, the shape and orientation of the halo should evolve on roughly a dynamical timescale. 
Because dynamical time decreases with density, the outer region of a halo will have signifi- 
cantly longer dynamical times than the central region, where the galaxy resides. Thus, the 
inner region (and associated galaxy) will acquire a misalignment from the tidal field more 
quickly than the halo as a whole. The process of misalignment has been studied in simula- 
tions by [38] , who find that halos orbiting in a cluster environment can experience significant 
internal twisting leading to appreciable misalignment between dark matter and a central 
stellar component. 

This connection of matter accretion to galaxy orientation yields two predictions. First, 
objects that have more recently accreted matter will be more strongly aligned with the tidal 
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field (corresponding to a larger value of Ci in the LA model). This trend is consistent with 
the observed luminosity dependence of the lA signal. Second, because of the gradient in 
dynamical time, stochastic evolution after the shut-off of the alignment mechanism, via grav- 
itational relaxation or interactions with other systems, will result in misalignment between 
the central regions of halos (or equivalently the galaxies they contain) and the outer halo 
regions, which will remain more closely oriented with the tidal field. Such misalignment is 
observed in simulations [28] with an amplitude consistent with that required to match align- 
ment statistics calculated from simulations with those observed from real galaxies [29, 30]. 
Using a semi-analytic galaxy catalog approach, [28] find that more luminous galaxies (as- 
sumed to have the same orientation as the central halo component) are more closely aligned 
with their parent halos. However, the trend they measure is too weak to account for all of the 
luminosity dependence in the alignment correlation function. Thus, halos containing more 
luminous objects may also exhibit stronger alignment with the tidal field, and galaxies may 
be less aligned with the central halo region than assumed. 

It would be interesting to more closely examine possible mechanisms for establishing 
intrinsic alignment. The interplay between accretion and other dynamical evolution may 
impart a redshift dependence on the strength of alignment. While accretion was generally 
more rapid at higher redshift, at later times the effects of additional accretion could increase 
alignment, or other evolution processes could become significant. Accurate characterization 
of this redshift dependence would improve our understanding of the importance of different 
physical processes in galaxy formation and evolution (although, see [32] who find no signif- 
icant redshift dependence for LRGs out to 2: ~ 0.7). Further measurements and theoretical 
work in this area are of particular relevance since the sources used in lensing studies are 
often at higher redshifts than those used to measure lA. Future work with both A^-body and 
baryon simulations can probe the dynamical aspects of galaxy alignment. Detailed observa- 
tions could provide a direct probe of misalignment between galaxies and the halos in which 
they reside. 

Encouraging progress has been made in recent years towards understanding the intrinsic 
alignment of galaxies and quantifying its effects. Much work remains to be done as we prepare 
for the next generation of weak lensing experiments. The issues discussed here, combined 
with the improved precision of future measurements, will provide methods to better describe 
and account for the intrinsic alignment of galaxies. 
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A Calculating w{rp,6) in the linear alignment model 

Consider a set of N Gaussian fields {qi}, which are assembled into a vector q, and some 
function of these fields F{qi, ...,q]y). An ensemble average of F is given by the integral: 

where Cij = (qiqj) is the relevant covariance matrix. 

We consider a reference galaxy at position with shear 70 and a second galaxy at point 
3^2 • The second galaxy has position angle 9, as measured from the major axis of the reference 
galaxy. Measured from the separation axis, the reference galaxy has 7+ = 70 cos (20). We are 
interested in calculating the 3-dimensional alignment correlation function, ig{r,9). We start 
with the definition of the standard correlation function, £,g{r): 

{nin2) = n^[l + Cg{r)], (A.2) 

where rii = n[l + 61] and from the 0-dependence, n2 = n[l + 62(0)]. Thus, 

{nin2\e) = n2[l + (Jal^) + (^i^al^)]. (A.3) 

Comparing with (A.2), we have: 

(g{r,9) = {62\e) + {5i52\e). (A.4) 

We apply eq. (A.l) with F = 62 + S162 and four Gaussian fields: {81,62, j+,^x}, 
where 7(+,x) is evaluated for the reference galaxy at x-^. Isotropy gives (7+7+) = (7x7x); 
(5i7(+_x)) = 0, and parity conservation requires (7+7x) = {52^x) = 0. These symme- 
tries greatly simplify the covariance matrix. Finally, we switch integration variables using 
d'y^d^jx = djodO det{ J) , where J is the Jacobian and det(J) = 270. Integrating over all 
variables except for 6, the result is 

\ 1/2 

^g{r, e) = ^g{r) + ( ^^J^^J cos(20)(527+)r, (A.5) 



where we have multiplied by an overall factor of vr to match the normalization convention 
adopted in the paper. Performing the projection yields 

Wg{rp,e) = Wg{rp) + ( 2(7_^_7+) ) cos{2e)wg+{rp). (A.6) 

We thus see that Gaussian density and ellipticity fields allow only cos (20) angular dependence 
in Wg{rp,6). From eqs. (3.9) and (A.6), we also find the following result: 

^ = (A.7) 

Wg+ 



TT 



where we have used 2(7_|_7+) = 7rjns. This ratio is the effective average of 70, where the 
ellipticity of each galaxy is weighted by how strongly it contributes to the statistic. If there 
were no correlation between 70 and 6, the ratio would simply be (70) for the entire galaxy 
sample. It is this hypothesis that is tested in [30], where they detect a ~ 15% modulation 
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effect due to a correlation between 70 and 9. Note that the results derived here (eqs. (A.5)- 
(A.7)) hold for any scenario in which the density and ellipticity fields are Gaussian, such as 
the LA model. 

We can employ the same formalism to more closely examine this correlation in the LA 
model. If we don't perform the 70 integral, we find: 

We can define an ellipticity- weighted alignment correlation function Wg{rp,6). Without the 
7o integral, this statistic is Wg{rp,6,'yo) = 70 w^g(?"p, 0, 7o)- 

The isotropic term gives the standard correlation function Wg{rp), while the angle- 
dependent term is related to Wg+. These terms have different 70 scalings because 82 and 
6162 pick out different terms in the exponent of the likelihood integral, eq. (A.l). The ratios 
between unweighted and weighted statistics will thus be different for isotropic and angle- 
dependent terms. Performing the integral over 70 yields 

Wg+{rp) TTWg{rp) ^ 



vr 



When we only consider counts of galaxy pairs, as with Wg{rp), all galaxies contribute 
equally. When pair counts are multiplied by alignment angle, as with galaxies that 

are more closely aligned with separation vectors are more heavily weighted. A correlation 
between 70 and 6 leads to a different effective average value for 79. For the LA model, such 
a correlation is present: more closely aligned galaxies have larger average values of 70. 

B Bessel function identities for E and B modes 

Inserting the LA model results into eq. (3.7), we find: 

'^^(^) =~2 (^ ^Diz) ^ j dKdk^-^^Ps{k,z)sm{kzIl^a.^) (B.l) 



Ja{^x) +4 r dx'^^i^ - 12^2 r dx 

J X J X 



x'^ 



Applying the recursion relation J„(y) = ^ [J„„i(y) + J„+i(y)] yields: 

My) = My) + - (-My) - My)) ■ (B.2) 
y \y J 

We then use following integral identities: 

Ja(kx') 

dx' \^ ' = [2J2(kj;) + HixMi^x)] , (B.3) 

dx' 

J X 

and find 

M^^x) + A r dx'^^^^ - Ylx" r dx'^^^^ = Jo(kx). (B.4) 

Jx Jx 

Thus, w'{x) = w+{x) in the LA model. 



, J^^KX') 


= 1 [ 


x' 




f J^^KX') 


_ Mi^x) 


x'^ 


KX^ 
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C Including redshift-space distortions 



When measured redshifts are used to determine distances to galaxies, peculiar velocities 
introduce errors in the line-of-sight position [55]. For the LA model, we are concerned with 
large scales, where coherent infall causes an apparent squashing of structure along the line- 
of-sight, leading to a spurious enhancement of clustering power. As first described in [55], 
we can relate the power spectrum in real space, Pr{k), with one measured in redshift space, 
Ps{k): 

Ps{k)=Pr{k)[l+(3fi^]^ (C.l) 

where /i = cos{9k) for angle 9^ between k and the line-of-sight; f3 = f{a)/bg for galaxy bias 
bg and f{a) = dlnD/dlna, the logarithmic growth rate of fluctuations. For this work, we 
apply the fitting formula /(a) = r2m(a)°'^ and values for bg measured using the amplitude of 
the correlation functions of the different galaxy samples. We use /3 ~ 0.3 for the LRG sample 
and /3 ~ 0.4 for the main sample, consistent with previous measurements [56]. Although we 
consider projected statistics, the effects of redshift-space distortions remain non-negligible 
because of the finite projection length (see, e.g., [57]). 
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